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We revisit the idea of using varying couplings to probe the nature of dark energy, in particular 
by reconstructing its equation of state. We show that for the class of models studied this method 
can be far superior to the standard methods (using type la supernovae or weak lensing). We also 
show that the simultaneous use of measurements of the fine-structure constant a and the electron- 
to-proton mass ratio yi allows a direct probe of grand unification scenarios. We present forecasts for 
the sensitivity of this method, both for the near future and for the next generation of spectrographs 
— for the latter we focus on the planned CODEX instrument for ESO's Extremely Large Telescope 
(formerly known as OWL) . A high-accuracy reconstruction of the equation of state may be possible 
all the way up to redshift z ~ 4. 



I. INTRODUCTION 

Observations suggest that the recent universe is dom- 
inated by an energy component whose gravitational be- 
haviour is quite similar to that of a cosmological constant 
(as first introduced by Einstein) 0, 0, 0] ■ This could of 
course turn out to be the cause for the recent acceleration 
of the Universe, but the observationally required value is 
much smaller than what would be expected from particle 
physics. In addition to the magnitude of the vacuum en- 
ergy, we are also faced with the cosmic coincidence that 
the vacuum energy density is of the same order of mag- 
nitude as the matter density. A dynamical scalar field 
is arguably a plausible explanation of the latter problem 
. Theoretical motivation for such a field is not hard to 
find. In string theory, for example, dimensional parame- 
ters are expressed in terms of the string mass scale and 
a scalar field vacuum expectation value. 

A crucial observational goal is therefore to characterize 
the properties of dark energy, and in particular to look 
for evidence of dynamical behaviour. A simple but im- 
portant property is its equation of state, w = p/p, and 
considerable effort has recently been put into trying to 
measure it as a function of redshift. The current meth- 
ods of choice are type la supernovae and (more recently) 
weak lensing, which is slightly more promising than the 
former. However, the question arises as to whether these 
are indeed the best tools for the task at hand. It has been 
known for some time [E ||J that supernova measurements 
are limited as a probe of the dark energy equation of 
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state, especially if it is varying with redshift. A discus- 
sion of current and future constraints on the dark energy 
equation of state, from the various standard approaches 
and parametrized in the usual way, is given in Ref. 0, 
and it is easy to conclude from it that a convincing de- 
tection of time variation of w is unlikely even with hypo- 
thetical future space-based experiments such as DUNE 
or JDEM (be it in its SNAP or DESTINY incarnations). 
This is not surprising since any dynamical field providing 
the dark energy must be slow-rolling at the present time 
(this is mandatory in order to have acceleration) , and for 
slow variations there will always be a constant w model 
that produces nearly identical results over the redshift 
range where dark energy is dynamically important. 

So, it is important to ask whether better (and cheaper) 
alternatives are available. A potentially effective tool for 
probing dynamical dark energy has been suggested pre- 
viously in |g, |2l Il0| , though not yet studied in detail. In 
any realistic model where the dark energy is due to a 
dynamical scalar field, that same field is also expected 
to produce sizable varying couplings ^E EE E3 ■ While 
some phenomenological studies of these models have been 

carried out Jjl EEM E3, [H, E3, H3, H3, H3, [S3, 1^, 

HE HE HE l2flj|. we show that probing these couplings can 
be a key test of these class of models, and in particular 
the varying couplings can be used to infer the evolution of 
the scalar field, and thus determine its equation of state. 
This is analogous to reconstructing the ID potential for 
the classical motion of a particle once its trajectory has 
been specified. 

Previous efforts [E IE EE only considered the variation 
of the fine-structure constant a, for which there is now 
a considerable dataset of observations and some (albeit 
disputed) evidence for variations in the redshift range 
z ~ 1 — 3 y(J, |3J, y^, |33| — see also |3J, y5| for reviews. 
However, in theories where a dynamical scalar field is 
responsible for varying a, the other gauge and Yukawa 
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couplings are also expected to vary |36|. Specifically, in 
GUTs there is a relation between the variation of a and 
that of the QCD scale, Aqcd, implying that the nucleon 
mass will vary when measured in units of the Planck 
mass. Similarly, one would expect variations in the Higgs 
vacuum expectation value (VEV), v, leading to changes 
in all particle mass scales including the electron mass. We 
therefore expect variations of the proton-to-electron mass 
ratio, \i = m^Jm^. Some specific models are discussed 
in Refs. [M El IS HI S3 S3 • Typically the relation 
between the variations is 
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the latter equality should be seen as the first term in 
a Taylor series, but given the expected level of varia- 
tions the approximation should be good enough for most 
purposes. The value of R is model-dependent (indeed, 
even its sign is not determined a priori), but large val- 
ues and negative values are generic for GUT models in 
which modifications come from high-energy scales: typi- 
cally A/A ~ 30a/a and v/v ~ 80d/a so that R ~ —50. 
The large proportionality factors arise simply because the 
strong coupling constant and the Higgs VEV run (expo- 
nentially) faster than a. 

If this intuitive picture is correct, then variations of // 
may be easier to detect than those of a. As first pointed 
out in Ref. > astrophysical measurements of fx can be 
made using Hi molecular lines. Some recent measure- 
ments using this technique 0, [13 a ls° nn d evidence for 
a variation of /i, although the number of current measure- 
ments is currently much smaller than the a dataset — the 
main reason for this is the difficulty in finding molecular 
Hydrogen clouds. One of the goals of the present paper 
is to encourage further measurements of /i, which are ex- 
pected to lead to tighter constraints on the evolution of 
the dark energy equation of state than those of a if R is 
indeed large. 



II. RECONSTRUCTING THE EQUATION OF 
STATE 

As we already pointed out, it was claimed that one 
could in principle extract information on the evolution of 
the equation of state of dark energy through the redshift 
dependence of a measurements @H, El • Here we briefly 
review the reconstruction procedure and show that in- 
cluding measurements of p in this approach allows us to 
estimate the value of the parameter R. We will concen- 
trate here on minimally coupled models of dark energy. 

The reconstruction pipeline follows the procedure in- 
troduced in Ref. for a. We assume that the functional 
form of the variation of these constants can be simply 
parametrized by a polynomial such that 
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where the various g xi are the coefficients of the poly- 
nomial, N = — ln(l + z) and x stands for either a or 
u. Given the tight Weak Equivalence Principle bounds 
|l4i Hsj , we also assume that the evolution of a and \x 
depends linearly on the value of the quintessence field 
through 
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and that £ M = R( a in agreement with Eq. (JTJ- 

Under these assumptions it was shown in Ref that 
the equation of motion for the energy density of the scalar 
field can be written as 
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where a = p^/po^Mo and a prime denotes differentiation 
with respect to N. The equation of state is obtained from 
the relation w = —1 + d> 2 1 p& which can be rewritten as: 
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Thus the evolution of the equation of state is determined 
by fitting the polynomial to the data, solving the 
differential equation (jJJ and by using its result in (j^J. 

The reconstruction under this procedure is ambiguous 
regarding the overall scale since the value of Q x is un- 
known. Hence a normalization of the equation of state 
must be performed by fixing its value at a given redshift 
to a value compatible with the value measured via in- 
dependent observations. More specifically, knowing the 
present values wq and fi^Oj we can estimate C, x through 
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which is nothing more than Eq. JSJ) at present (noting 
that 1 + l/aa 3 — l/O^). Once this coupling has been 
estimated we can proceed by solving Eq. therefore 
obtaining the evolution of the equation of state at large 
redshifts only using information based on the variation 
of constants. In practice current observations only con- 
strain effective values of fi^ and w (a weighted average 
over redshift) rather that their present day values. How- 
ever, our simpler assumption will not greatly affect our 
results. We also stress that for any given ( x it is possible 
to obtain any evolution for Ax/x by choosing appropriate 
scalar field potential and initial conditions. This means 
that (a; can only be effectively constrained using comple- 
mentary datasets. 

If the same dynamical scalar field providing the dark 
energy is also responsible for the varying couplings there 
will be an important consistency test provided by viola- 
tions of the Equivalence principle. Indeed, the smaller 
Q x is, the larger (j) needs to be in order to be responsi- 
ble for a given evolution of Ax/x. Hence, Q x cannot be 
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made arbitrarily small or otherwise the kinetic energy of 
the field would be too large for it to be the dark energy. 
In this context, it has been shown [23l liH lirj that, if 
at least some of the claimed detections of varying cou- 
plings are correct then there must be violations of the 
Einstein Equivalence Principle at a level within reach of 
the forthcoming generation of space-based experiments 
(ACES, AiSCOPE, GG and STEP), which is expected to 
improve on the present sensitivity by as much as five or- 
ders of magnitude. Hence, one can envisage that these 
experiments might themselves soon provide a measure- 
ment of £ x rather than just bounds. 

There is also the possibility of combining our recon- 
struction with a standard method such as using type la 
Supernovae. Despite its weaknesses, which we discussed 
above, reconstruction using type la Supernovae is a more 
direct probe of the evolution of 6 and wihprovide an im- 
portant consistency check [H El El Efl El El In 
particular, we should be able to test, though with lim- 
ited accuracy, the validity of our main assumption (the 
linear relation between Ax/x and </> given by Eq. ©) 
using future space-based experiments such as DUNE or 
JDEM. 



III. SIMULATED DATA 

We perform an analysis having in mind the quality of 
the datasets expected to be available both in the near 
future (say in 2 or 3 years time) and what is expected 
for the next generation of ground-based astronomical in- 
strumentation. Specifically, for the latter case we focus 
on the expected sensitivity of the CODEX instrument for 
the European ELT (formerly known as OWL). 

Note that we will restrict our attention to direct mea- 
surements of the various couplings. The analysis of the 
Oklo natural nuclear reactor (see 0] for a critical re- 
view) offers very accurate measurements of a at z ~ 0.14, 
Aa/a ~ 1CT 7 . Further more this limit is strengthened by 
about two orders of magnitude if one assumes that other 
couplings are varying |4£j. Nevertheless, these methods 
rely on data (and theory) beyond our direct control and 
as such will not be applied here. Another example of 
an indirect measurement is one based on OH lines in 
gravitational lens systems at a redshift of order z ~ 0.7 
[54j. On the other hand, measurements using the CMB 

[H, EI El El 

are direct but share with Oklo and also 
with local laboratory experiments E3| the disadvan- 
tage of measuring a mainly at a single redshift, which 
makes the inter-comparisons subtle (due to the differing 
systematics), as well as model-dependent. Methods that 
can probe a wide range of redshifts are therefore best 
suited for our purposes. 

We produced simulated data sets for both a and /i 
using a Monte Carlo generator based on the redshift de- 
pendence of the scalar field for three particular scalar 
potentials under a corresponding choice of the parame- 
ter £ x . The evolution of the field with redshift was de- 



termined by numerically integrating the field's equation 
of motion for the chosen potentials. These are models 
previously investigated in the literature and we fix their 
corresponding parameters so that the various examples 
can be representative of the possible evolutions of the 
equation of state. Two of the models are given by a 
potential with two exponential terms [6J. Model A, is 
given by V(<p) = Vo(exp(lOK0) + exp(0.1/c<^>)) and model 
B by V(<f>) = Uo(exp(50K</>)+exp(0.8K(/>)). The first term 
allows the energy density of the field to approach an 
attractor solution for a wide range of initial conditions 
and the second term drives the Universe into a period 
of accelerated expansion in accordance with current ob- 
servational constraints. In the first model the equation 
of state parameter w, increases with redshift and in the 
second it decreases. The third model, model C, is given 
by, V{4>) = V Q exp(2(n(t>) 2 ) |2j|. For large values of the 
scalar field, this potential also has an attractor solution 
with w ~ constant, however, the equation of state pa- 
rameter quickly decreases to —1 as the field approaches 
the minimum. 

Based on the current data and the theoretical esti- 
mates for the parameter R, we consider a fiducial model 
in which the parameter Q. is fixed is such a manner that 
we obtain (Aa/a) 2=3 = —0.06 x 10~ 5 and (A/i/ fj) z =3 = 
3 x 10~ 5 , thus R = —50. For models A, B and C, it can 
be found that the appropriate values are £ Q = 1.6 x 10~ 6 , 
( a = 2.07 x 10~ 6 and ( a = -0.74 x 10~ 6 , respectively. 

The distributions of the data sets are equivalent to 
the current Murphy et al. data [3l| . This data spans 
the interval z = [0.2,3.7], however, we interpolate the 
distributions up to redshift z — 4.2 as a damped Lyman- 
a (DLA) absorption system has been recently identified 
at this redshift. 

We assume that in the near future we will be able to 
determine the variation in a from a sample of 250 objects 
and the variation in \x from 20 objects. We take the min- 
imum value of the error bars to be w 10~ 6 occurring at 
a redshift z » 1.5 and increasing quadratically from this 
redshift. The coefficient of this quadratic function and 
the redshift of minimum error bar were also determined 
from the present data set. An example of the typical 
simulated data sets is shown in the upper two panels of 
Fig. □ For the CODEX spectrograph the precision m 
a is estimated to be « 10~ 8 . We use this value for the 
minimum error bar in the simulated data. As before, 
the error bars of the remaining data follow a quadratic 
function with minimum near redshift z = 1.5. We conser- 
vatively use a sample of 500 data points for a and of 100 
data points for \i. The high resolution simulated data is 
shown in the lower two panels of Fig. ^ 



IV. RESULTS 

We have performed a likelihood analysis by fitting 
the various parameters g xi to the simulated data sam- 
ples. When using both samples together, we fit g^- 
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FIG. 1: Potential data sets, generated from model A, ex- 
pected in the relatively near future (top two panels labeled 
with S > 10~ 6 ) and with CODEX (lower two panels labeled 
with S > 10~ 8 ) as described in the text are compared. S rep- 
resents the minimum value of the error bars considered in the 
corresponding simulated data. 



and in addition R such that the full vector sample is 
[R(Aa/a)j a , (A/i//i) Jt J, where j a and are indices run- 
ning over the size of the a and /x samples, respectively. 



A. Near-future prospects 

In Figs. [3 El an d [J] w e show the potential for recon- 
structing the equation of state parameter, w as a function 
of redshift. The dark (light) shaded regions correspond 
to the 1 (2) a confidence levels of reconstruction. We 
clearly observe from the upper left panels of Figs. [21 El 
and 0] that the reconstruction based entirely on Act /a 
will still be unsatisfactory in the near future. The recon- 
struction shown in the upper right panel of Figs. El Eland 
0]of the equation of state based on the fi sample, however, 
will be very illuminating. For the fiducial model under 
study, 20 points are quite sufficient to make it clear that 
w increases with redshift. Given that the signal in fi is 
fifty times larger than the one on alpha, it is no surprise 
that the reconstruction is in this case more successful. 
We can read from Tables □ |TT| and EH that we can ob- 
tain an improvement on the reduced \ 2 going to second 
order in the degree of the polynomial g^, though a third 
parameter is only helpful when reconstructing model B. 
Because the data on a is comparatively so poor, when 
combining both samples we obviously find that the pa- 
rameter R is also poorly constrained (see Fig. EJ). 

Alternatively, one could consider generating a simula- 
tion based on the result from the analysis in Ref. |3l| 
implying that (Aa/a) z= 3 w —0.5 x 10~ 5 (in which case 
we would have R ~ —5). In this case, we would conclude 
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TABLE I: Best fit results when fitting the a samples and 
fi samples with m g xi parameters, and both samples with 
the additional parameter R. The left columns shows the 
reduced x 2 an d Ax 2 = X 2 ( m — 1) — X 2 ( m )j f° r the type 
of data expected in the near future and the right columns 
the same quantities for data expected with CODEX. The 
data was generated from the evolution using the potential 
V((j>) = Vo(exp(10«</>) +exp(O.lfi:0)) (model A). 
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TABLE II: Same as in Table □ for the potential V((f>) 
Vo(exp(50K</>) +exp(O.8«0)) (model B). 
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TABLE III: Same as in Tables Q] and E] for the potential 
V(<t>) = Voexp(2(K^) 2 ) (model C). 
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FIG. 2: The reconstruction to the mth order of the equation 
of state and its error band is shown for the various surveys 
assuming the potential V(<f>) = Vb(exp(10ft</>) + exp(0.1n(j>)) 
(model A). The dashed line represents the dark energy equa- 
tion of state corresponding to the potential used to generate 
the simulated data and the solid line corresponds to the re- 
construction's best fit. The dark region is the la confidence 
level on the parameters and the light region the 2a confidence 
level. 
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FIG. 3: Same as in Fig. with the potential V{4>) = 
Vb(exp(5O«;0) +exp(O.8«0)) (model B). 



that the a sample, though not as good as a /i sample, 
provides nonetheless a comparable reconstruction with a 
similar error band and the accuracy of the determination 
of R can be significantly better than the one presented 
in Figs. El El and E| 
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FIG. 4: Same as in Figs. Eland El with the potential V(4>) = 
V exp(2(K(t>) 2 ) (model C). 



B. The promise of CODEX 

We now turn our attention to the reconstruction of 
the equation of state derived from measurements with 
the CODEX instrument In this case, the sam- 

ples allow for a precise reconstruction and provide the 
possibility to go beyond the first coefficients of the poly- 
nomial. From Tables HI ITU and InTl we see that a fit for a 
involving at least three parameters is now preferable. In 
principle the error bars are small enough to admit a poly- 
nomial of even higher degree (see lower panels of Figs. El 
El and El This will be in fact necessary if the quintessence 
field oscillates at the minimum of its potential. There 
are, however, difficulties associated with models with os- 
cillatory behaviour as explained in Ref. [9| . These models 
have typically an equation of state near — 1 which intro- 
duces a large uncertainty in the value of Q x and conse- 
quently on the overall reconstruction of the equation of 
state. However, as discussed above this difficulty might 
be avoided with future datasets. Combining both data 
sets, the uncertainty on R will be limited to only a frac- 
tion of a percent as can be shown in Fig. El 

Consider a sample with N a and independent data 
points for Aa/a and Afj,/ fj, with a given redshift distribu- 
tion. Let us also assume a given redshift distribution for 
the individual Aa/a and Afx/n error bars with a nor- 
malization given respectively by e Q and e M . The error 
bands in the reconstruction of the equation of state scale 
approximately as 



where N a and correspond to the size of the a and /i 
samples respectively. In fact Eq. Q allows for a simple 
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FIG. 5: The precision achieved on the parameter R with the 
various surveys by fitting a and [t, samples together generated 
from model A. The solid line illustrates the possible likelihood 
to be obtained with near future data and the dashed line, the 
likelihood using the CODEX data. 

re-scaling of the results if different samples with simi- 
lar redshift distributions for the Aa/a and Afi/fi data 
points and respective error bars are used. For the sim- 
ulated data used to perform the reconstruction of the 
equation of state our simple re-scaling formula gives a 
ratio of 14:1:0.1:0.004 for the relative sizes of the error 
bands of the various models considered in Figs. [3 El and 
2] which is roughly in agreement with the above results. 

V. DISCUSSION 

We have discussed the possible use of varying couplings 
to reconstruct the equation of state of dark energy, and 
emphasized the high benefit of obtaining further mea- 
surements of fi. In combination with measurements of 
a (which are easier to obtain, given the abundance of 
sources) they can provide a key test of the relationship 
between varying couplings in a grand unification sce- 
nario. This method not only complements results antici- 
pated by hypothetical future experiments (such as JDEM 
and DUNE), but given reasonable expectations for forth- 
coming improvements in spectroscopic measurements, is 
expected to be competitive with the standard methods 
for dark energy equation of state reconstruction (both 
those using supernovae and those based on weak lens- 
ing). 

Having two different observables (the fine-structure 
constant and the proton-to-electron mass ratio) , one can 
in fact check the self-consistency of the reconstruction. In 
fact, given that standard procedures directly probe the 
evolution of <j> at low redshift, another consistency check 



will be provided by combining our equation of state re- 
construction with the standard one. A failure of the con- 
sistency check can be interpreted in one of two ways: that 
the coupling @ is inadequate and possibly an expansion 
to higher order in cj) is in order; or that the coupling to 
matter vanishes altogether and it is impossible to corre- 
late the variations in couplings with dark energy. 

For the class of models studied, however, the recon- 
struction of the dark energy equation of state using vary- 
ing couplings has several fundamental advantages over 
the standard methods. Firstly, one has the advantage 
of a much larger lever arm in terms of redshift, since 
such measurements can easily be made up to redshifts of 
z ~ 4. This may not seem like a substantial advantage, 
since dark energy is only dynamically relevant at rela- 
tively low redshift, but in fact it is a key one, since one 
can probe the redshift range where the field evolution is 
expected to be fastest (if it is a tracking field) — that is, 
deep in the matter era. Secondly, the varying couplings 
method is also much cheaper and can be done without 
any problems from the ground, even with existing facili- 
ties - all it takes is a few hundred good nights of telescope 
time, which is quite a modest investment given the po- 
tential gains. Note that the optimal observing strategy 
(e.g., the choice of the number of sources to observe and 
their redshift distribution), for a given total available ob- 
servation time, will depend on a range of factors, both 
theoretical and observational. The design of such a strat- 
egy is an important issue that needs to be studied in more 
detail in the future. 

To conclude, let us again stress the power and simplic- 
ity of this method. It requires only ground-based obser- 
vations, and a relatively small dataset (possibly requir- 
ing not more than a hundred hours of telescope time) 
can have a huge impact and conceivably provide unam- 
biguous evidence for dynamical dark energy. The main 
point to keep in mind is that a good and uniform red- 
shift coverage is important, since one will have to cal- 
culate derivatives of the processed data. The small but 
representative range of models we have considered shows 
that this is within reach in the next few years. Last 
but not least, this is also an example of how astrophys- 
ical observations can be optimal probes of fundamental 
physics. We believe that such astrophysical probes will 
become increasingly common in years to come, and hope 
this provides early encouragement for the observational 
astrophysics community. 
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